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Density functional theory based global geometry optimization has been used to demonstrate the crucial influence of the geometry 
of the catalytic cluster on the energy barriers for the CO oxidation reaction over Pd-based bimetallic nanoalloys. We show that 
dramatic geometry change between the reaction intermediates can lead to very high energy barriers and thus be prohibitive for 
the whole process. This introduces challenges for both the design of new catalysts, and theoretical methods employed. On the 
theory side, a careful choice of geometric configurations of all reaction intermediates is crucial for an adequate description of 
a possible reaction path. From the point of view of the catalyst design, the cluster geometry can be controlled by adjusting 
the level of interaction between the cluster and the dopant metal, as well as between the adsorbate molecules and the catalyst 
cluster by mixing different metals in a single nanoalloy particle. We show that substitution of a Pd atom in the Pds cluster with 
a single Ag atom to form Pd 4 Agi leads to a potential improvement of the catalytic properties of the cluster for the CO oxidation 
reaction. On the other hand, a single Au atom does not enhance the properties of the catalyst, which is attributed to a weaker 
hybridization between the cluster’s constituent metals and the adsorbate molecules. Such flexibility of properties of bimetallic 
nanoalloy clusters illustrates the possibility of fine-tuning, which might be used for design of novel efficient catalytic materials. 


1 Introduction 

Understanding catalysis at the atomic level is one of the fun¬ 
damental goals of chemistry. Supported metal cluster catalysis 
in particular plays an enormous role in the modern chemical 
industry, and much effort is devoted to understanding this phe¬ 
nomenon and to the design and development of metal-cluster- 
based heterogeneous catalysts with high activity, selectivity, 
and stability. 1 Bimetallic catalysts often exhibit advantageous 
properties compared to those of their pure constituent metals, 
as this hybrid material may offer opportunities for the syn¬ 
ergistic inter-metallic interactions that lead to improvement 
in catalytic performance. 2 Moreover, the catalytic activity of 
bimetallic alloys or nanoparticles can be optimized by control¬ 
ling structural factors, such as the alloying element and con¬ 
centration. 3 The geometric and electronic effects have been 
shown to systematically alter the catalytic activity of bimetal¬ 
lic catalysts. 4 

For example, it has been found that the addition of metals 
that bind oxygen strongly (e.g. Co, Ni, and Cu) to more noble 
metals (Pt or Pd) can improve their oxygen reduction catalytic 
activity. 5-8 Bimetallic clusters are also known to be effective 
catalysts for methane dehydrogenation, 9 and a rich variety of 
organic reactions, such as the Suzuki reaction, 10 the dehalo- 
genation of aromatic cycles, 11,12 and the oxidation of benzyl 
alcohol, 13 among others. 
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A very important process catalysed by bimetallic clus¬ 
ters is the CO oxidation reaction, which has a large impact 
in terms of environmental protection. Great efforts have 
been devoted to developing highly efficient catalysts to re¬ 
move the poisonous CO gas from car exhaust and from hy¬ 
drogen gases used for fuel cells. 14 Theoretical studies of 
the CO adsorption and oxidation have been carried out for 
Au-Pd, 15 Pd-Pt, Cu-Pt, Pd-Rh, 16 Zr-Sc, 17 and Au-Cu 18 
clusters. Experimentally, catalytic CO oxidation over Pd-Ni 
has been reported recently. 19 Theoretical works propose both 
Langmuir-Hinshelwood (two molecules adsorb and the ad¬ 
sorbed molecules undergo a bimolecular reaction) and Eley- 
Rideal (one of the molecules adsorbs and the other one reacts 
with it directly from the gas phase) mechanisms. For example, 
in the case of small Au-Pd clusters, both reaction mechanisms 
were investigated, with the Langmuir-Hinshelwood found to 
be more energetically favourable. 15 

If there is oxygen present in the system, such as in the case 
of Zr-Sc oxide clusters, oxygen atoms for CO oxidation are 
provided by the cluster. 17 In the case of Au-Cu clusters ad¬ 
sorbed on the oxide surface, the co-adsorption of CO and O 2 
(Langmuir-Hinshelwood) takes place. 18 Thus, the Langmuir- 
Hinshelwood mechanism seems to be preferable for CO oxi¬ 
dation on bimetallic clusters. Another concern is the stability 
of the catalyst upon adsorption. The data in Ref. 16 indicates 
the possibility of structural rearrangements in the Au-Pd clus¬ 
ters upon adsorption of the CO. 

Theoretical methods are therefore especially suited for 
studying the catalytic properties of such complex materials, 
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as they allow one to obtain detailed information on the mech¬ 
anism of the reaction and the structural rearrangements in¬ 
volved, which is not directly accessible to experiment. More¬ 
over, theory can aid the design of a novel catalyst by testing 
the catalytic activity of the systems of variable size and com¬ 
position against a suitable performance descriptor . 20 

A number of theoretical studies have been conducted to elu¬ 
cidate the influence of the dopant metals on the catalytic activ¬ 
ity of clusters . 21,22 It has been suggested that substitution of a 
single Pd atom in a Pds cluster with a Au atom might improve 
the catalytic activity . 21 Moreover, it has been suggested that 
adsorption of such a Pd4Au cluster on a TiO2(110) surface 
leads to a more efficient CO to CO2 conversion, compared 
to a pure Pds cluster. At the same time, pure gold clusters 
have also been recently proposed to potentially exhibit high 
catalytic activity for the same process, e.g. Aug adsorbed on 
ZnO . 22 Thus the question of the optimal ratio of Pd and Au 
atoms in the catalytic Pd-Au clusters remains unsolved. 

A very important open problem is finding the most rele¬ 
vant configurations of the intermediates at every stage of the 
reaction. The reaction path obviously depends on the cluster 
configuration, as concominantly do the heights of the energy 
barriers. Ideally one would like to characterize the complete 
potential energy surfaces for clusters and adsorbates to iden¬ 
tify the most likely pathway for the catalytic reaction. How¬ 
ever, due to the expense of the level of theory needed to obtain 
reliable results and the large number of different structures one 
has to consider for alloy clusters, this is potentially not feasi¬ 
ble. This is the main bottleneck of a theoretical investigation 
of catalysis on nanoalloys. What is the best approach in the 
absence of a complete characterization of a potential energy 
surface is an open question. 

Often the catalyst is assumed stable during the reaction, and 
only small adsorbate molecules (such as O2 and CO in our 
case) are allowed to undergo major configurational changes, 
while the geometry of the catalytic metal cluster is only re¬ 
laxed locally. This assumes that the rearrangement of the clus¬ 
ter geometry is slow compared to the overall rate of catalytic 
reaction. Assuming such stability of the cluster is likely to 
lead to relatively simple reaction paths and perhaps relatively 
low transition state energy barriers. However, if there is suffi¬ 
cient energy and time to allow significant reconfigurations in 
the adsorbate geometry, why should the cluster geometry not 
relax to more stable geometries if available? 

An alternative approach is to assume that the cluster adopts 
the ground-state configuration at every step in the reaction. 
This approach assumes that the rearrangement of the cluster 
is relatively fast compared to the overall rate of the catalytic 
reaction, and requires a global optimization of the geometry 
at every step. To validate the applicability of this approach, 
such a “ground-state-based” reaction path should, of course, 
be compared to other possible paths involving low-lying iso¬ 


mers of every intermediate. 

These open questions motivate the present theoretical study 
of the CO oxidation catalysed by Pd-based bimetallic nanoal¬ 
loys. We first find a reasonable Pd/dopant ratio based on the 
adsorption energies of CO and O2 molecules on the corre¬ 
sponding clusters. After that we study the complete reaction 
path of the CO oxidation for the Pd clusters alloyed with all 
possible dopants from the typically catalytically active groups 
VIIIB and IB: Au, Ag, Cu, Ni, Pt. Here, we rely on the ther¬ 
modynamic feasibility of the catalytic process, thus effectively 
using the heights of the energy barriers along the reaction path 
as a descriptor of whether the catalyst is potentially active. Fi¬ 
nally, we discuss the observed differences in the behaviour of 
the clusters of different composition based on an analysis of 
the electronic structure. 


2 Computational details 

All local geometry optimizations of the discussed structures, 
and subsequent electronic structure analysis were carried out 
with the plane-wave DFT package CASTEP. 23 Electronic ex¬ 
change and correlation was treated within the generalized- 
gradient approximation functional due to Perdew, Burke and 
Ernzerhof (PBE ). 24 The core electrons were described us¬ 
ing ultrasoft pseudopotentials, whereas the valence electrons 
were treated with a plane-wave basis set with a cut-off en¬ 
ergy of 300 eV. Local structure optimization is done us¬ 
ing the Broyden-Fletcher-Goldfarb-Shanno method , 25 relax¬ 
ing all force components to smaller than 10 -2 eV/A. The sta¬ 
bility of the identified minima has been confirmed by vibra¬ 
tional frequency analysis. 

To obtain the ground-state structures for all considered 
systems, we relied on basin-hopping (BH) based global ge¬ 
ometry optimization , 26,27 using the DFT total energies and 
atomic forces calculated by CASTEP 23 as implemented in the 
Atomic Simulation Environment (ASE) suite . 28 

The Nudged Elastic Band (NEB ) 29 method as implemented 
in the ASE 28 package was used to find the transition paths and 
corresponding energy barriers between given initial and final 
states, which were defined as the ground-state structures of 
adjacent reaction intermediates. Permutational isomers of the 
global minima were chosen in such a way as to minimize the 
geometry change involved in the transition. A chain of twelve 
replicas of the system has been constructed for every transition 
state calculation, and the forces were allowed to relax with a 
threshold of 0.05 eV. The two states are assumed to be con¬ 
nected by a single barrier path. For the transitions involving 
the loss of CO2 the corresponding final minimum is reopti¬ 
mized with a desorbed CO2 molecule explicitly present in the 
system. 
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The reaction path of CO oxidation over a bimetallic catalyst 
is defined as 

2 CO -f- O2 “l - PdyA/y —y 2CO2 “l - P d x M y5 

where M is a doping metal atom. The possible inter¬ 
mediate products are assumed as CO-Pd x M v , 02-Pd x My, 
0 2 C 0 -Pd x My, O-Pd x M y , and OCO-Pd x M } , The total 
exothermic effect of the reaction is computed to be 6.41 eV 
( 147.8 kcal/mol), with the experimental value reported as 
135 kcal/mol. 30 

3 Selecting suitable cluster size and Pd/dopant 
ratio 

In order to identify a system that is likely to have good cat¬ 
alytic properties, we suggest that a candidate catalyst should 
meet the following requirements: 

• Adsorption of CO and O2 is reasonably strong and of 
similar strength. 

• Dissociative adsorption of O2, or a substantial destabi¬ 
lization of the 0—0 bond is observed. 

• Reactions of the adsorbed species on the cluster are char¬ 
acterized by low barriers. 

• All the intermediates are less stable than the products of 
the reaction. 

• The cluster shows little affinity for CO2. 

The first issue to address is to find a suitable size of the clus¬ 
ter, and a suitable ratio of atom types in the cluster. We use the 
first two requirements from the above list to assess different 
clusters. Ref. 21 suggested a catalytically active Pd5 cluster 
as a reference point for Pd to Au substitution. Following this 
logic, we decided to either substitute atoms in such a cluster, 
thus ending up with a set of six Pd5_ x Au x clusters, or to cap it 
with an additional atom, thus having a set of seven Pd6- x Au x 
clusters of various composition. We started with Pd-Au clus¬ 
ters as a model system, and later introduced other transition 
metal dopants for nanoalloys with the selected Pd/dopant ratio 
that was thought most likely to yield good catalytic properties. 

As can be seen in Fig. 1 , clusters having more than two 
gold atoms have a tendency to prefer planar geometries as 
their ground states. Thus, AU5, PdiAu4, Au6, and PdiAus 
are planar both as individual clusters, and with either CO or 
O2 adsorbed. This is in line with the known tendency of small 
gold clusters to retain planar configurations up to the size of at 
least AU7, 31 and at even larger sizes in the case of anionic clus¬ 
ters. 32 Pd 2 Au 3 and Pd2Au4 clusters turn into planar structures 
upon adsorption of CO or O2. This can be explained in terms 
of the reduced hybridization between the states of Pd and Au, 
when the Pd atoms predominantly interact with an adsorbate 
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Fig. 1 Lowest-energy Pd5_ x Au x and Pd 6 _ x Au x clusters, and their 
aggregates with CO and O2. Numbers indicate adsorption energies 
in eV. 


molecule rather than with the neighbouring Au atoms. Indeed, 
rather strong interactions between the adsorbates and the Pd 
atoms of the cluster (up to 2.49 eV binding energy in the case 
of the CO-Pd2 Au 3 aggregate) allow the relatively undisturbed 
gold atoms to retain the relativistic aspects 33 of their prop¬ 
erties. This is also reflected in the tendency of gold atoms 
to cluster together in these cases. Another interesting feature 
worth mentioning is the fact that both CO and O2 tend to di¬ 
rectly interact with Pd atoms, and not Au atoms. The Pd 3 Au 3 
cluster can be qualified as an intermediate case, showing a pla¬ 
nar configuration for the CO-Pd 3 Au 3 aggregate, and a com¬ 
pact geometry for the 02-Pd 3 Au 3 aggregate. 

Clusters with a fraction of gold atoms that is less than 
half possess compact structures, and, more importantly, have 
higher binding energies towards both CO and O2. In cer¬ 
tain cases (such as Pd2Au 3 and PdiAu4) the O2 molecule is 
not even dissociated upon adsorption. Thus a sufficient frac¬ 
tion of Pd atoms allows more efficient activation of the ad¬ 
sorbate molecules for the oxidation reaction, while compact 
structures ensure the stability of the cluster itself under realis¬ 
tic conditions, such as deposition on a support surface. From 
analysis of the CO and O2 adsorption energies we conclude 
that Pd4 Aui and Pd4 AU2 represent particularly good composi¬ 
tions in terms of the adsorption efficiency: these cluster sizes 
have compact geometries, ensure the dissociation of the O2 
molecule, and have the highest binding energies towards CO 
and O2 among all considered clusters. 
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This conclusion is also backed up by similar calculations 
performed for the clusters deposited on the TiO2(110) sur¬ 
face, which is a widely used support for catalytic metal clus¬ 
ters in the reaction of CO oxidation. 34-36 Pd4Aui and Pd4Au2 
provide the highest adsorption energies for both CO and O2 
molecules compared to other Pd/Au ratios, except for pristine 
Pd clusters: deposited Pds and Pd6 yield the highest adsorp¬ 
tion energies in all cases. Fig. 2 depicts the ground-state struc¬ 
tures for the pure clusters and their aggregates with CO and 
O2 in the cases of Pd4Aui and Pd4Au2. For comparison, the 
corresponding structures of pristine Pds and Pds clusters are 
also provided. 



Fig. 2 Lowest-energy Pd 5 , Pd 4 Aui, Pds, and PdsAui clusters 
adsorbed on the Ti02 (110) surface, and their aggregates with CO 
and O 2 . Only the top layer of the Ti02 support is displayed for 
clarity. Numbers indicate adsorption energies in eV. See 
supplemental material 37 for the structures and energies of other 
Pd/Au compositions. 

Interestingly, the oxygen atoms in the top layer of the 
TiO2(110) surface can participate in the stabilization of the 
clusters with adsorbed molecules, as can be seen by the change 
in their positions compared to the clean surface. By con¬ 
trast, bare clusters retain their stability upon deposition, and 
do not disturb the surface much. However, in the presence of 
the adsorbate molecules the picture is quite different. While 
the ground-state geometries of the CO-Pd5 and CO-Pd4Aui 


aggregates are very similar to their pure cluster counterparts, 
the addition of the CO adsorbate molecule changes the global 
minima configurations of the deposited Pd6 and Pd4Au2 clus¬ 
ters significantly. The adsorption of the O2 molecule leads to 
rather distorted final structures in all cases. Even more striking 
is the fact that the addition of gold atoms prevents dissociation 
of the oxygen molecule. For both Pd4Aui and Pd4Au2 clus¬ 
ters, the ground-state structures of the aggregates with oxygen 
correspond to activation of the O2 molecule (0-0 bond length 
of 1.32 A and 1.35 A in the case of Pd4Aui and Pd4Au2, re¬ 
spectively, compared to 1.21 A for a gas-phase O2 molecule), 
but not dissociation. This behaviour can be explained by the 
observed partial charge transfer from the deposited cluster to¬ 
wards the surface. The resulting electron deficit prevents the 
cluster from becoming an active site for the dissociation of 
the adsorbed O2 molecule which has a high electron affin¬ 
ity. However, previous theoretical studies 21,38-40 suggest that 
this fact does not necessarily prevent the TiC>2 ( 110 )-deposited 
clusters from exhibiting efficient catalytic properties. 

4 Selecting suitable dopants 

Thus in our further investigations we focus on the cluster of 
and Pd4M2. As dopant atoms we have chosen the met¬ 
als from the typically catalytic groups VIIIB and IB, i.e. Au, 
Ag, Cu, Ni, and Pt. For comparison, pristine Pd5 and Pd6 clus¬ 
ters were also considered. In order to investigate the detailed 
mechanism of the reaction we initially take the “global mini¬ 
mum approach” where we assume all intermediates adopt their 
lowest-energy structure. Consequently, we run global geom¬ 
etry optimization for every intermediate of the assumed reac¬ 
tion path, i.e. Pd x M y , CO-Pd x M y , 02-Pd x M y , 02C0-Pd x M y , 
O-Pd x M y , and OCO-Pd x M y . Transition state structures were 
found by applying the NEB method using the ground-state 
structures of the adjacent intermediates as initial and final con¬ 
figurations, respectively. Thus the first transition state (TS 1 ) 
corresponds to the reaction path between the C^CO-Pd x M y 
and 0 -Pd x My with a first CO2 molecule released. The sec¬ 
ond transition state (TS 2 ) corresponds to the reaction path be¬ 
tween OCO-Pd x My and an individual cluster with a second 
CO2 molecule released. Note that the final configuration of 
the cluster (after two CO2 molecules have been produced) is 
assumed to match the initial ground-state structure of the clus¬ 
ter. In the following we first consider the reaction paths for 
Pds-based and Pd6-based clusters separately, and then discuss 
the common features. 

4.1 Reaction paths for Pds-based clusters 

The identified reaction paths for the Pds-based clusters are 
presented in Fig. 3 . The reaction path for the pristine Pds 
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Fig. 3 Reaction paths for Pds-based clusters: (a) Pd 5 , (b) Pd 4 Aui, (c) Pd 4 Agi, (d) Pd 4 Cui, (e) Pd 4 Nii, (f) Pd 4 Pti. 


cluster (Fig. 3 (a)) already exhibits the main features of a rea¬ 
sonably good catalyst. Firstly, every further intermediate has 
a relative energy which is lower than that of a previous one. 
Secondly, the transition state (TS) energy barriers are rela¬ 
tively low ( 1.38 eV for TS 1 and 0.96 eV for TS 2 , see Table 1 ). 
Furthermore, both CO-Pds and O2-PCI5 aggregates have very 
similar relative energies, which makes it more likely that both 
CO and O2 bind to the cluster with similar probability, thus 


preventing the cluster from being predominantly covered by 
the more strongly binding species. Also noteworthy is that 
the global geometry optimization revealed that the geometry 
of the cluster itself does not change much during the reac¬ 
tion. In such cases, the barriers reflect just those structural 
changes needed for the adsorbates to react, and not any addi¬ 
tional structural change of the cluster. 

Substitution of one Pd atom with a Au atom (Fig. 3 (b)) 
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does not change the energy profile of the reaction drastically. 
However, there is an important difference in how the clus¬ 
ter geometry changes during the reaction. In particular, the 
lowest-energy structures of the intermediates involving simul¬ 
taneous adsorption of the CO and either an O2 molecule or an 
O atom (i.e. 02C0-Pd4Aui at — 4.42 eV and OCO-Pd4Aui 
at — 6.07 eV in Fig. 3 (b)) have a quasi-planar geometry. Such 
geometries turn out to be energetically most favourable, and 
yield lower relative energies, compared to the corresponding 
intermediates in the case of the pristine Pds cluster. There¬ 
fore a substantial rearrangement of the cluster atoms is conse¬ 
quently needed to get to the next reaction intermediate, which 
probably underlies the higher energy barriers compared to Pd5 
(Table 1 ). As a result, the overall efficiency of the catalytic CO 
oxidation is (somewhat surprisingly) expected to be lower in 
the case of a gold-substituted Pd4Aui compared to the pure 
Pds cluster. 

Substitution with a silver atom instead of a gold atom 
(Fig. 3 (c)) helps to overcome the problem of the cluster ge¬ 
ometry changing during the reaction. As can be seen from 
the reaction path for the Pd4Agi cluster, none of the inter¬ 
mediates exhibit planar cluster geometries, thus the reaction 
involves less structural rearrangements and has lower energy 
barriers. Especially striking is the difference for the last step 
of the reaction OCCMAUMi —» Vd/^M\ + CO2, where in the 
case of Pd4Agi, only the adsorbate atoms undergo significant 
rearrangement; this step involves a very low energetic barrier 
of 0.55 eV, compared to the 1.36 eV for Pd4Aui. Thus, the 
Pd4Agi cluster appears to be a better potential catalyst for 
CO oxidation than Pd4Aui. Experimentally, Pd-Ag bimetal¬ 
lic clusters have been recently shown to be catalytically active 
for various reactions. 41 ^ 3 

Doping with a Cu atom also leads to reasonably good cat¬ 
alytic properties (Fig. 3 (d)). With copper ensuring compact 
structures for all the intermediates, somewhat similar to those 
of their Ag-substituted counterparts, the reaction path of CO 
oxidation over the Pd4Cui cluster features a second energy 
barrier ( 1.09 eV for TS 2 ) that is lower than that of the Pd4Aui 
cluster. However, the first energy barrier ( 1.65 eV for TS 1 ) 
is comparable to that of a gold-substituted cluster. Although 
no substantial structural change is involved, the relatively high 
energy barrier can probably be attributed to the higher affinity 
of oxygen atoms towards copper compared to silver. It is also 
worth mentioning that CO2 is able to physisorb on the Pd4Cui 
cluster (note that for consistency the last energy in the reac¬ 
tion profile in Fig. 3 (d) corresponds to the total energy release 
of the CO oxidation reaction, assuming CO2 dissociated from 
the cluster). However, the small adsorption energy of 0.16 eV 
should not lead to any noticeable catalyst poisoning at finite 
temperature under actual experimental conditions. Good per¬ 
formance of Pd-Cu nanoalloys (although for a much larger 
size of -5 nm) in the process of CO oxidation has also been 


Table 1 Transition state energy barriers for Pd 5 -based clsuters 


Cluster 

TS1 

Energy barrier, eV 

TS2 

Pd 5 

1.38 

0.96 

Pd4Aui 

1.57 

1.36 

Pd 4 Agi 

1.34 

0.55 

Pd4Cui 

1.65 

1.09 

Pd 4 Nii 

1.86 

1.69 

Pd 4 Pti 

1.83 

2.48 


observed experimentally. 44 

Using Pt or Ni as dopants for a Pd cluster does not lead to a 
potentially effective catalyst. In both cases (Fig. 3 (e) and (f)) 
the last OCO-Pd4Mi intermediate is more stable than the cor¬ 
responding metal cluster with a newly formed CO2 molecule, 
thus preventing the release of the second CO2, and leading to 
poisoning of the catalyst. This tendency is also apparent from 
the clusters with O2 and CO adsorbate molecules: both clus¬ 
ters have significantly higher stability for the C^CO-PcUAfi 
intermediate than either the pristine Pds cluster, or clusters 
with Au, Ag, and Cu dopant atoms. This effect is most pro¬ 
nounced in the case of Pd4Nii cluster: here the affinity for O2 
is so high that the difference between the adsorption energies 
of CO and O2 is as large as 1.45 eV. The 02C0-Pd4Nii in¬ 
termediate is also more stable than the final product of the 
CO oxidation reaction (— 6.47 eV vs. — 6.41 eV), and the 
02C0-Pd4Ni -> 0-Pd4Ni + CO2 reaction leading to the re¬ 
lease of the first CO2 is unfavourable by 0.91 eV. This can be 
explained by the higher mobility of the Ni ^-electrons, where 
the 3 d^ 4 s 2 configuration allows active donation towards ox¬ 
idizing agents, thus preventing Ni from exhibiting properties 
similar to its group VIIIB noble metal neighbours. 

4.2 Reaction paths for Pd6-based clusters 

Overall, Pd6-based clusters exhibit similar CO oxidation be¬ 
haviour. However, there are important differences caused by 
the increased complexity of the systems. For example, the 
rather strong structural distortion, which the pure Pd6 cluster 
undergoes in TS 1 , changing from a 02C0-Pd6 intermediate 
to a 0 -Pd 6 aggregate (Fig. 4 (a)), is reflected in a rather high 
energy barrier of 2.10 eV (Table 2 ). The second transition 
(OCO-Pd6 Pd6 + CO2) does not involve much structural 
change, and yields a smaller 1.35 eV barrier. 

Substitution of two Pd atoms with two Au atoms en¬ 
hances the catalytic activity compared to the pure Pd6 clus¬ 
ter (Fig. 4 (b)). Unlike in the case of the singly-doped Pd4Aui 
cluster, all the Pd4 Au2 intermediates are non-planar. Although 
somewhat counterintuitive, the increased fraction of gold in 
the doubly-doped cluster compared to its singly-doped coun¬ 
terpart does not lead to an increased probability of forming 
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CO oxidation over a Pd 4 di 2 cluster 





Fig. 4 Reaction paths for Pd6-based clusters: (a) Pd6, (b) Pd 4 Au 2 , (c) Pd 4 Ag 2 , (d) Pd 4 Cu 2 , (e) Pd 4 Ni 2 , (f) Pd 4 Pt 2 . 


planar ground-state structures for the intermediates. How¬ 
ever, a noticeable rearrangement of the cluster’s geometry is 
required during the first transition, which has a rather high 
energy barrier of 1.78 eV, compared to 1.57 eV for the cor¬ 
responding transition in the singly-doped Pd4Aui cluster. By 
contrast, the second transition exhibits a much smaller distor¬ 
tion and yields a somewhat lower energy barrier ( 1.29 eV for 
Pd 4 Au 2 vs. 1.36 eV for Pd4Aui). As a result, Pd4Au2 overall 


shows more favourable properties as a potential catalyst than 

Pd 6 . 

Silver dopants further increase the efficiency of the doubly- 
doped Pd 4 Ag 2 cluster catalyst (Fig. 4 (c)). Smaller structural 
distortions are reflected in the decreased values of transition 
state energy barriers of 1.58 eV for TS 1 , and 1.17 eV for TS 2 . 
Although the properties of Pd4 Ag2 are thus potentially better 
than those of either Pd4Au2 or Pd6, the Pd4Ag2 has barriers 
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Table 2 Transition state energy barriers for Pd^-based clsuters 


Cluster 

TS1 

Energy barrier, eV 

TS2 

pd 6 

2.10 

1.35 

Pd 4 Au 2 

1.78 

1.29 

Pd 4 Ag 2 

1.58 

1.17 

Pd 4 Cu 2 

2.01 

1.50 

Pd 4 Ni 2 

1.50 

2.68 

Pd 4 Pt 2 

2.98 

2.94 


that are larger than both its smaller Pd4Aui counterpart, and 
the Pds cluster. 

The Cu-substituted cluster (Fig. 4 (d)) exhibits energy bar¬ 
riers similar to those of the pure Pd6 cluster. However, the 
last OCO-Pd4Cu2 intermediate is 0.56 eV more stable than a 
recovered catalyst cluster with a second CO2 molecule pro¬ 
duced. Moreover, physisorption of CO2 that was noted earlier 
for the singly-copper-doped cluster is even more pronounced 
here: the 0.38 eV adsorption energy for Pd4Cu2 is a sign of 
chemical bond formation. 

As is also the case for clusters, Pt and Ni dopants 

do not assist the formation of a potentially effective Pd4M2 
catalysts (Fig. 4 (e) and (f)). Addition of a second Ni atom 
enhances the affinity towards oxygen atoms even further: the 
difference between the adsorption energies of O2 and CO is 
as large as 2.05 eV. This affinity to oxygen is also reflected in 
the value of the second transition state barrier: 2.68 eV ob¬ 
viously makes the reaction unfeasible. In the case of Pd4Pt2, 
additional Pt atoms cause very strong distortions to the cluster 
geometry, which is reflected in very high energy barriers of 
almost 3 eV for both transition states. 

It should be pointed out that that many of the higher en¬ 
ergy barriers are associated with transitions between interme¬ 
diates with distinctively different cluster structures. This can 
be indicative of the difference in chemical nature of Pd clus¬ 
ters doped with different metals, or the increased complex¬ 
ity of the system, e.g. Pd4M2 compared to Pd4Mi. The key 
question, however, is whether there are alternative pathways 
that lead to an overall lower reaction barrier, e.g. a multi-step 
process where a cluster might first change its geometry, and 
then the reaction occur. Such pathways may be missed by our 
single-path NEB calculations. 

We suggest that such methodological limitations will not 
change our conclusion that the Ni- and Pt-doped clusters (and 
also Pd 4 Cu 2 ) are unlikely to be good catalysts, since the un¬ 
favourable poisoning of these clusters with CO2 is indepenent 
of whether there are alternative lower-energy paths for the re¬ 
actions. For the other nanoalloys it is desirable to consider 
different possible reaction paths. This problem is addressed in 
section 4 . 4 . 


4.3 Electronic structure analysis of Au- and Ag-doped 
clusters 

Comparison of the CO oxidation efficiency on Pd4Mi and 
Pd 4 M 2 clusters reveals that doping can be an effective way 
of tuning the catalytic properties of Pd clusters. Thus, dop¬ 
ing with silver lowers the energy barriers in both cases. The 
behaviour of gold dopants is, however, more complex. While 
Pd 4 Au 2 exhibits potentially better catalytic properties than its 
pristine Pd6 counterpart, the smaller Pd4Aui yields higher 
transition state energy barriers due to the observed tendency 
to form planar intermediate structures, which leads to signifi¬ 
cant geometric rearrangement during the reaction steps. 

In order to explain the observed differences between Au and 
Ag dopants, and also to look into the nature of chemical bond¬ 
ing of both CO and O2 towards catalytic clusters, we calcu¬ 
lated density of states diagrams for the Pd4Aui and Pd4Agi 
cluster. As can be seen in Fig. 5 , stronger bonding between 
the dopant and Pd, as well as strong binding of the adsorbates 
is observed in the case of the silver dopant. Kohn-Sham levels 
of silver and palladium show pronounced hybridization with 
each other in the region between — 10.0 and — 7.5 eV (Fig. 
5 (b)). Moreover, a strong contribution from the adsorbates is 
also observed in this region. While a smaller interaction is also 
observed between palladium and gold (Fig. 5 (a)), no strong 
hybridization with the adsorbates is observed. 

Insets in Fig. 5 depict the highest occupied molecular or¬ 
bitals (HOMO) of Pd4Aui and Pd4Agi, which allows the be¬ 
haviour of the most chemically relevant Kohn-Sham states 
to be directly examined. Since both clusters are triplets in 
their ground state, the HOMOs correspond to the highest oc¬ 
cupied orbitals in the alpha spin channel. The HOMO of 
Pd4Aui lies at — 5.39 eV with respect to the vacuum level, 
while the HOMO of Pd4Agi is situated at — 5.16 eV. The dif¬ 
ferences between the dopants are best observed in the case of 
CO dopant adsorption. As was also observed in a recent the¬ 
oretical study, 21 the bonding between the CO molecule and 
Pd4Aui is mostly of the s-p cr-type, with a large fraction of 
electron density on the CO molecule. In the case of a sil¬ 
ver atom dopant, however, a more complicated ^-character is 
observed, indicating stronger interaction of the CO molecule 
with the cluster. The C—O bond length in the case of Pd4Agi 
is increased to 1.19 A, compared to 1.16 A in Pd4Aui, and 
1.14 A in equilibrium), which suggests that CO is better acti¬ 
vated by the silver-doped cluster. 

Overall, weaker adsorbate-cluster interaction, and a smaller 
degree of metal-metal hybridization for an Au dopant indi¬ 
cates that doping with Au reduces the reactivity of the cluster, 
compared to the addition of Ag. This renders silver as a poten¬ 
tially better dopant for the design of novel Pd-based bimetallic 
catalysts for CO oxidation. 
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Fig. 5 Electronic densities of states projected on Pd (black line), the dopant atom (red line), and adsorbate molecules (blue line) for (a) 
02C0-Pd4Aui, and (b) 02C0-Pd4Agi. The vacuum level is used as the zero reference. The insets show the highest occupied molecular 
orbital (HOMO) of the corresponding cluster. 


4.4 Relative TS barriers for reaction vs. rearrangement 

Calculation of the reaction paths based on the globally opti¬ 
mal structures for each intermediate assumes that adopting the 
ground-state structure of each intermediate is kinetically more 
likely than the immediate reaction of the adsorbate molecules 
on the surface of an unperturbed cluster. To explore whether 
this assumption is reasonable, we have chosen two systems for 
which pronounced changes in cluster geometry were observed 
along the reaction pathway, namely Pd4Aui and Pd4Pt2. For 
each intermediate we then analyse low-lying isomers identi¬ 
fied during the global optimization run, and choose the one 
corresponding to the least perturbation of the cluster geome¬ 
try relative to the ground-state structure of the individual clus¬ 
ter. We then construct the reaction path diagram for such an 
“alternative” or “unperturbed” reaction path, and compare it 
to the “ground-state-based” reaction path diagram created on 
the basis of the globally optimal intermediates. Additionally, 
we also analyse the energy barriers for the reconfiguration of 
each reaction intermediate structure from the “unperturbed” to 
its ground-state geometry in order to understand whether such 
reconfiguration is likely to take place during the reaction. 

As was suggested by the analysis of the ground-state-based 
reaction path, reconfigurations of the Pd4Aui reaction inter¬ 
mediates are reflected in a slight increase in the transition state 
energy barriers for such singly-doped cluster compared to the 
pristine Pd5 cluster. To investigate whether the “alternative” 
reaction path leads to lower energy barriers, in Fig. 6 we plot 
both the energy levels of the “alternative” and the ground- 
state-based reaction paths. As can be seen from the Fig. 6, 
all “unperturbed” intermediates are naturally less stable than 


CO oxidation over a Pd 4 Aui cluster 



Fig. 6 Possible reaction paths of the CO oxidation over Pd 4 Aui: 
“alternative” reaction path (black lines) compared to the 
ground-state-based reaction path (grey lines). Insets correspond to 
the “unperturbed” intermediates. 


the ground-state structures. Thus, at the first step of the reac¬ 
tion, the CO-Pd4Aui aggregate is 0.36 eV less stable than the 
global minimum, while 02-Pd4Aui lies 0.62 eV higher than 
its globally optimal counterpart. However, the reconfiguration 
of these structures into the global minima is likely: the cor¬ 
responding downhill energy barriers are only 0.06 eV in both 
cases (Fig. 7 ). 

The energetic difference between the “unperturbed” and 
the ground-state configurations of the C^CCMAUAui is 
0.63 eV. However, reconfiguration between these quite dif- 
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ferent structures is energetically expensive: a 2.25 eV bar¬ 
rier needs to be overcome. Although the first reaction barrier 
(02C0-Pd4Aui 0-Pd4Aui + CO2) is slightly lower in the 
case of the “alternative” reaction path ( 1.38 eV vs. 1.57 eV), 
the transition state involved is 0.44 eV higher in energy than 
the transition state for the ground-state-based path, and so it is 
still likely to be more favourable to proceed along the ground- 
state-based path. 

The next intermediate, 0-Pd4Aui, lies 0.23 eV higher in 
energy than its ground-state structure with a downhill energy 
barrier of a mere 0.03 eV, which makes the formation of the 
global minimum configuration very likely. The two isomers 
of the final intermediate, OCO-Pd4Aui, lie very close to each 
other in energy: there is only 0.17 eV difference between the 
two. Although the transition to the global minimum requires a 
rather high barrier of 1.52 eV to be surmounted, which makes 
the transition less favourable. The “alternative” OCO-Pd4Aui 
isomer has a higher energy barrier for the CO oxidation reac¬ 
tion compared to the ground-state-based path, namely 1.54 eV 
compared to 1.36 eV. 

For the CO-Pd4Aui, 02-Pd4Aui, and 0-Pd4Aui interme¬ 
diates the reconfiguration required for a transition to the corre¬ 
sponding global minima is a simple diamond-square-diamond 
process between trigonal pyramidal isomers with the Au atom 
in apical and equatorial positions. This process facilitates 
an easy transition towards the more energetically favourable 
ground-state-based reaction path at the early stage of the re¬ 
action. By contrast, in the case of the 02C0-Pd4Aui and 
OCO-Pd4Aui intermediates, the transition between an “un¬ 
perturbed” cluster and the corresponding global minimum in¬ 
volves transformations between a non-planar geometry and a 
planar one, which is energetically costly, and which may make 
the escape from the initially chosen pathway unlikely. 

In the case of Pd4Pt2, the ground-state-based reaction path 
features very high energy barriers, and leads to the poison¬ 
ing of the catalyst due to high affinity of Pd-Pt clusters to¬ 
wards oxygen. Can the “alternative” reaction path based on 
the unperturbed pure cluster geometry be more energetically 
favourable? As can be seen from the Fig. 8, perturbation of 
the structure is a natural result of the interaction of the nanoal¬ 
loy cluster with the adsorbate molecules, especially the O2 
molecule: both first reaction intermediates, i.e. CO-Pd4Pt2 
and O2—Pd4Pt2 are 0.63 eV and 1.60 eV less stable than their 
ground-state structures. 

Such a significant energetic gain from following the 
ground-state-based reaction path is also conserved in the case 
of the second reaction intermediate C>2CO-Pd4Pt2 ( 1.63 eV 
less stable than the global minimum) and the first transi¬ 
tion state TS1 (02C0-Pd4Pt2 -A 0-Pd4Pt2 + CO2), which 
lies 2.19 eV higher in energy than its ground-state-based path 
counterpart. This yields a very high energy barrier of 3.54 eV. 
The resulting 0-Pd4Pt2 intermediate is in turn significantly 
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Fig. 7 Pd 4 Aui intermediate aggregates reconfiguration barriers. 


less stable than its previously identified ground-state isomer 
( 1.41 eV higher in energy). 

Upon adsorption of the second CO molecule, 0-Pd4Pt2 
turns into OCO-Pd4Pt2, which in this case is only slightly less 
stable than the ground-state structure ( 0.35 eV difference). 
However, this energy is already lower than the total effect of 
the reaction by 0.28 eV, which again means poisoning of the 
catalyst, as was also observed in the case of the ground-state- 
based reaction path. 

The second reaction transition state energy barrier turns 
out to be significantly lower for the “alternative” re¬ 
action path ( 1.07 eV vs 2.94 eV), as only adsorbate 
molecules undergo reconfiguration during the last reaction 
step (OCO-Pd4Pt2 —» Pd4Pt2 + CO2). However, both the pro¬ 
hibitively high first transition state energy barrier, and the ob¬ 
served poisoning of the catalyst are common to both reaction 
pathways and do not allow the success of the catalytic pro¬ 
cess. The energy barriers associated with the reconfiguration 
between the intermediates of the ground-state-based and the 
alternative reaction paths are presented in the supplemental 
material. 37 

The above results support the initially proposed significant 
role of the geometric configurations of the reaction interme¬ 
diates on the observed transition state barriers. In the case 
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CO oxidation over a Pd 4 Pt 2 cluster 



Fig. 8 Possible reaction paths of the CO oxidation over Pd 4 Pt 2 : 
“alternative” reaction path (black lines) compared to the 
ground-state-based reaction path (grey lines). Insets correspond to 
the “unperturbed” intermediates. 

of Pd4Aui, the ground-state-based reaction path is a reason¬ 
able assumption, as the corresponding intermediate structures 
are more stable, and the transition between the higher lying 
“unperturbed” isomers and their ground-state counterparts are 
kinetically likely for the first two reaction steps. In the case 
ofPd 4 Pt 2 , much less stable intermediate geometries still lead 
to a significant distortion of the structure at the stage of the 
release of the first C 0 2 molecule, thus yielding an increased 
transition state energy barrier. Moreover, the final reaction 
intermediate is still lower in energy than the reaction prod¬ 
ucts, which corresponds to the poisoning of the catalyst. Thus 
choosing ground-state configurations of the reaction interme¬ 
diates appear to yield energetically reasonable reaction paths 
for the considered systems. 

Notwithstanding, the example of Pd4Pt 2 illustrates the com¬ 
plexity of the configurational space that needs to be sampled 
to arrive at a definitive conclusion, and hence the potential pit- 
falls of our theoretical approach. As already mentioned above, 
these include the danger of picking sub-optimal configurations 
for the reaction intermediates, and a potential tendency of the 
employed transition state search method to miss more compli¬ 
cated multi-step pathways connecting two intermediates that 
might potentially involve smaller barriers. 

5 Conclusions 

In summary, we have systematically studied the influence of 
different transition metal dopants on the catalytic activity of 
Pd-based clusters in the reaction of CO oxidation. Unbiased 
DFT-based basin-hopping global geometry optimization re¬ 
veals that the transition state energy barriers of the reaction 


are very sensitive to the geometric configurations of the in¬ 
termediate structures. For example, substituting a Pd atom in 
the Pds cluster with a Au atom leads to the increase of both 
transition barriers due to the tendency of gold-doped clus¬ 
ters to form planar structures. Noticeable rearrangement of 
the cluster geometry is thus needed to proceed along the reac¬ 
tion path, which is associated with substantial energetic costs. 
On the other hand, doping with silver facilitates the formation 
of compact structures, and helps to lower the energy barri¬ 
ers. Electronic structure analysis indicates that gold and silver 
dopants exhibit different chemical behaviour: silver atoms al¬ 
low stronger hybridization with the electronic levels of both 
the Pd atoms in Pd cluster and the adsorbate molecules, thus 
creating more favourable conditions for catalytic activation of 
such molecules. This renders silver a potentially better dopant 
for the CO oxidation reaction than previously suggested gold. 

Intriguingly, copper-doped nanoalloys show a potential for 
catalytic activity that is similar to that of the pristine palladium 
clusters. However, the higher affinity of oxygen for copper 
eventually leads to the poisoning of the catalyst in the case 
of Pd4Cu 2 . Bimetallic clusters with Ni and Pt do not exhibit 
favourable catalytic properties for CO oxidation due to high 
affinity for oxygen (Ni) and substantial structural rearrange¬ 
ments between the reaction steps (Pt). 

Size effects might also play an important role in such 
bimetallic systems. Our calculations indicate that the addi¬ 
tion of the second dopant atom increases the complexity of 
the system, and makes stronger geometric distortions between 
the reaction steps more likely, which is in turn linked to an 
increase in the transition state energy barriers. As a result, for 
both Au and Ag dopants the smaller Pd4Mi clusters appear 
to have better catalytic properties than the larger Pd4M 2 ones. 
From the methodological perspective, one has to take into ac¬ 
count that the increased configurational complexity also raises 
the cost of a thorough sampling of the potential energy surface, 
thus making it less likely that we have found the lowest-energy 
pathway and making it more important that additional calcu¬ 
lations considering alternative reaction paths are undertaken. 

Investigation of alternative reaction paths for the examples 
of Pd4Aui and Pd4Pt 2 clusters revealed that the ground-state- 
based reaction path, where each reaction intermediate is found 
as a result of the global geometry optimization, is likely to lead 
to a pathway with lower-energy transition states compared to 
the “unperturbed” reaction path, where each intermediate is 
represented by a low-lying isomer corresponding to the least 
perturbation of the cluster geometry compared to the ground 
state of the bare cluster. However, this is not guaranteed, par¬ 
ticularly if there are large changes in geometry between the 
global optima for the reactants and products of a catalytic step. 
In the absence of a rigorous study of all feasible reaction paths 
and their possible intersections, which is prohibitively expen¬ 
sive for the number of systems considered here, relying on 
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global geometry optimization for finding the configurations of 
the reaction intermediates appears to be a reasonable choice 
for constructing an energetically favourable reaction path. 

Overall, the study reflects the rich diversity of bimetal¬ 
lic cluster structures, and high tunability of their geometric 
and electronic properties with a choice of a transition metal 
dopant, as well as the cluster size and composition. This ren¬ 
ders such bimetallic nanoalloys as ideal candidates for the de¬ 
sign of novel catalytic materials. We would also like to em¬ 
phasize the methodological challenge of an adequate descrip¬ 
tion of the configurational space of bimetallic clusters, which 
makes a rigorous investigation of all possible reaction paths at 
the electronic structure theory level unfeasible. The task ide¬ 
ally calls for the development of general, reliable, and com¬ 
putationally efficient methods of describing the potential en¬ 
ergy surface of bimetallic clusters and adsorbate molecules, 
perhaps based on empirical potentials, but such a task is of 
course very challenging. 
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I. Structures and relative energies of the Pd5_ x Au x and Pd 6 - X Au^ clusters 

adsorbed on the TiC>2 (HO) surface 
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FIG. SI. Lowest-energy Pds-^Au^ and Pde-^Au^ clusters adsorbed on the Ti02 (110) surface, 
and their aggregates with CO and O 2 . Numbers indicate adsorption energies in eV. 
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II. Pd4Pt2 intermediate aggregates reconfiguration 
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FIG. S2. Pd 4 Pt 2 intermediate aggregates reconfiguration barriers. 









